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The problem of flame propagation is studied as an example 
of unstable fronts that wrinkle on many scales is studied. The 
analytic tool of pole expansion in the complex plane is em- 
loyed to address the interaction of the unstable growth process 
with random initial conditions and perturbations. We argue 
that the effect of random noise is immense and that it can 
never be neglected in sufficiently large systems. We present 
simulations that lead to scaling laws for the velocity and ac- 
celeration of the front as a function of the system size and 
the level of noise, and analytic arguments that explain these 
results in terms of the noisy pole dynamics. 
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I. INTRODUCTION 

The aim of this paper is to examine the role of ran- 
dom fluctuations on the dynamics of growing wrinkled 
interfaces which are governed by non-linear equations of 
motion. We are interested in those examples for which 
the growth of a flat or smooth interface are inherently un- 
stable. A famous example of such growth phenomena is 
provided by Laplacian growth patterns [1-3]. The exper- 
imental realization of such patterns is seen for example 
in Hele-Shaw cells [1] in which air or another low vis- 
cosity fluid is displacing oil or some other high viscosity 
fluid. Under normal conditions the advancing fronts do 
not remain flat; in channel geometries they form in time 
a stable finger whose width is determined by delicate ef- 
fects that arise from the existence of surface tension. In 
radial geometry, the growth the interface forms a con- 
torted and ramified fractal shape. A related phenomenon 
has been studied in a model equation for flame propa- 
gation which has the same linear stability properties as 
the Laplacian growth problem [4]. The physical problem 
in this case is that of pre-mixed flames which exist as 
self-sustaining fronts of exothermic chemical reactions in 
gaseous combustion. Experiments [5] on flame propaga- 
tion in radial geometry show that the flame front accel- 
erates as time goes on, and roughens with characteristic 
exponents. Both observations did not receive proper the- 
oretical explanations. It is notable that the channel and 
radial growth are markedly different; the former leads to 
a single giant cusp in the moving front, whereas the latter 
exhibits infinitely many cusps that appear in a complex 



hierarchy as the flame front develops [6,7]. 

Analytic techniques to study such processes are avail- 
able [8]. In the context of flame propagation [7,9-11], 
and in Laplacian growth in the zero surface-tension limit 
[12-14] one can examine solutions that are described in 
terms of poles in the complex plane. This description 
is very useful in providing a set of ordinary differential 
equations for the positions of the poles, from which one 
can deduce the geometry of the developing front in an 
extremely economical and efficient way. Unfortunately 
this description is not available in the case of Laplacian 
growth with surface tension, and this makes the flame 
propagation problem very attractive. However, it suf- 
fers from one fundamental drawback. For the noiseless 
equation the pole-dynamics always conserves the num- 
ber of poles that existed in the initial conditions. As a 
result there is a final degree of ramification that is af- 
forded by every set of initial conditions even in the ra- 
dial geometry, and it is not obvious how to describe the 
continuing self-similar growth that is seen in experimen- 
tal conditions or numerical simulations. Furthermore, as 
mentioned before, at least in the case of flame propaga- 
tion one observes [5] an acceleration of the flame front 
with time. Such a phenomenon is impossible when the 
number of poles is conserved. It is therefore tempting to 
conjecture that noise may have an important role in af- 
fecting the actual growth phenomena that are observed 
in such systems. In fact, the effect of noise on unsta- 
ble front dynamics has not been adequately addressed in 
the literature. From the point of view of analytic tech- 
niques noise can certainly generate new poles even if the 
initial conditions had a finite number of poles. The sub- 
ject of pole dynamics with the existence of random noise, 
and the interaction between random fluctuations and de- 
terministic front propagation are the main issues of this 
paper. 

We opt to study the example of flame propagation 
rather than Laplacian growth, simply because the for- 
mer has an analytic description in terms of poles also in 
the experimentally relevant case of finite viscosity. We 
choose to limit the present study to channel geometry. 
The reason is that in radial geometry it is more difficult 
to disentangle the effects of external noise from those 
of initial conditions. After all, initially the system can 
contain infinitely many poles, very far away near infinity 
in the complex plane (and therefore having an infinitely 
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small contribution to the interface). Since the growth 
of the radius changes the stability of the system, more 
and more of these poles might fall down to the real axis 
and become observable. In channel geometry the analysis 
of the effect of initial conditions is relatively straightfor- 
ward, and one can understand it before focusing on the 
(more interesting) effects of external noise [9]. The ba- 
sic reason for this is that in this geometry the noiseless 
steady state solution for the developed front is known 
analytically. As described in Section II, in a channel of 
width L the steady-state solution is given in terms of 
N(L) poles that are organized on a line parallel to the 
imaginary axis. It can be shown that for any number 
of poles in the initial conditions this is the only attrac- 
tor of the pole dynamics. After the establishment of this 
steady state we can begin to systematically examine the 
effects of external noise on this solution. As stated be- 
fore, in radial conditions there is no stable steady state 
with a finite number of poles, and the disentanglement of 
initial vs. external perturbations is less straightforward 
[7]. We believe nevertheless that the insights provided 
in this paper have relevance for radial growth as well as 
will be discussed in the final Section and in forthcoming 
publications. 

We have a number of goals in this paper. Firstly, after 
introducing the pole decomposition, the pole dynamics, 
and the basic steady state, we will present a stability 
analysis of the solutions of the flame propagation prob- 
lem in a channel geometry. It will be shown that the 
giant cusp solution is linearly stable, but non-linearly 
unstable. These results, which are described in Section 
III, can be obtained either by linearizing the dynamics 
around the giant cusp solutions in order to study the 
stability eigenvalues, or by examining perturbations in 
the form of poles in the complex plane. The main result 
of Section III is that there exists one Goldstone mode 
and two modes whose eigenvalues hit the real axis pe- 
riodically when the system size L increases. Thus the 
system is marginally stable at particular values of L, and 
it is always nonlinearly unstable, allowing finite size per- 
turbations to introduce new poles into the system. This 
insight allows us to understand the relation between the 
system size and the effects of noise. In Section IV we 
discuss the relaxation dynamics that ensues after start- 
ing the system with "small" initial data. We study the 
coarsening process that leads in time to the final solution 
of the giant cusp, and understand from this what are the 
typical time scales that exist in our dynamics. We offer 
in this Section some results of numerical simulations that 
are interpreted in the later sections. In Section V we fo- 
cus on the phenomenon of acceleration of the flame front 
and its relation to the existence of noise. In noiseless con- 
ditions the velocity of the flame front in a finite channel is 
bounded [9]. This can be shown either by using the pole 
dynamics or directly from the equation of motion. We 
will present the results of numerical simulations where 



the noise is controlled, and show how the velocity of the 
flame front is affected by the level of the noise and the 
system size. The main results are: (i) Noise is responsible 
for introducing new poles to the system; (ii) For low lev- 
els of noise the velocity of the flame front scales with the 
system size with a characteristic exponent; (iii) There is 
a phase transition at a sharp (but system-size dependent) 
value of the noise-level, after which the behavior of the 
system changes qualitatively; (iv) After the phase transi- 
tion the velocity of the flame front changes very rapidly 
with the noise level. In the last Section we remark on the 
implications of these observations for the scaling behav- 
ior of the radial growth problem, and present a summary 
and conclusions. 



II. EQUATIONS OF MOTION AND 
POLE-DECOMPOSITION IN THE CHANNEL 
GEOMETRY 

It is known that planar flames freely propagating 
through initially motionless homogeneous combustible 
mixtures are intrinsically unstable. It was reported that 
such flames develop characteristic structures which in- 
cludes cusps, and that under usual experimental condi- 
tions the flame front accelerates as time goes on. A model 
in 1 + 1 dimensions that pertains to the propagation of 
flame fronts in channels of width L was proposed in [4]. 
It is written in terms of position h(x, i) of the flame front 
above the «-axis. After appropriate rescalings it takes 
the form: 



dh(x, t) 
dt 



dh(x, t) 



1 2 



da 



d 2 h(x,t) 
dx 2 



+ I{h(x,t)} + 1 



(1) 



The domain is < x < L and v is a parameter. The func- 
tional I[h(x, t)] is the Hilbert transform which is conve- 
niently defined in terms of the spatial Fourier transform 



h(x, t) = 
I[h(k,t)} 



e lkx h{k,t)dk 

■ CO 

\k\h(k,t) 



(2) 
(3) 



For the purpose of introducing the pole-decomposition it 
is convenient to rescale the domain to < 6 < 2tt. Per- 
forming this rescaling and denoting the resulting quanti- 
ties with the same notation we have 
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dh(6,t) 
dt 

+jl{h(6,t)} + l 



dh(6,t) 
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v d 2 h(6,t) 
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dO 2 



(4) 



In this equation L = L/2tt. Next we change variables to 
u(6,t) = dh(6,t)/d6. We find 
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du(e,t) _ u(6,t) du(6,t) v d 2 u(6,t) 1 



dt 



L 2 
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L 2 



dO 2 



-I{u{9,t)} ■ 
(5) 



It is well known that the flat front solution of this equa- 
tion is linearly unstable. The linear spectrum in k- 
representation is 



ui k = \k\/L-vk 2 /L 2 



(6) 



There exists a typical scale k max which is the last unsta- 
ble mode 



L 



(7) 



Nonlinear effects stabilize a new steady-state which is 
discussed next. 

Following [8,9,11,7] we expand the solutions u(9,t) in 
functions that depend on N poles whose position Zj (t) = 
xj(t) + iyj(t) in the complex plane is time dependent: 
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u(6A 



) = "£ 



cot 



0- ZJ (t) 



+ c.c. 



N 



2sin[6» - Xj{t)] 



J cosh[t/j(t)] — cosf^ — xj(t)] 



(8) 



terms of the real and imaginary parts xj and yj . Because 
of the arrangement in pairs it is sufficient to write the 
equation for either yj > or for yj < 0. We opt for the 
first. The equations for the positions of the poles read 



~dt 



N 

v sm( Xj 



Xk) 



[cosh(t/ i - y k ) (11) 



■ cos(;Ej - x k )] 1 + [cosh(t/j + y k ) - cos(xj - x k ) 
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4- y 



sinhfa- - y k ) 



di k Jrt^j ^ cosll (%' ~ Vk) ~ cos(xj - x k ) 



sinh(t/ i + y k ) 



cosh(t/j + y k ) - cos(xj - x k ) 



) + i^coth(t/ i ) - L. (12) 



We note that if the initial conditions of the differential 
equation (5) are expandable in a finite number of poles, 
these equations of motion preserve this number as a func- 
tion of time. On the other hand, this may be an unstable 
situation for the partial differential equation, and noise 
can change the number of poles. This issue will be exam- 
ined at length in Section V. We turn now to a discussion 
of the steady state solution of the equations of the pole- 
dynamics. 



N 

h(9, t) = 2v In [ cosh( % - (t j) - cos(9 - Xj (*))] + C(t) . 
i=i 

(9) 

In (9) C(t) is a function of time. The function (9) is a 
superposition of quasi-cusps (i.e. cusps that are rounded 
at the tip). The real part of the pole position (i.e. xj) 
is the coordinate (in the domain [0, 2ir]) of the maximum 
of the quasi-cusp, and the imaginary part of the pole 
position (i.e yj) is related to the depth of the quasi-cusp. 
As yj decreases the depth of the cusp increases. As yj — ► 
the depth diverges to infinity. Conversely, when yj — ► 
oo the depth decreases to zero. 

The main advantage of this representation is that the 
propagation and wrinkling of the front can be described 
via the dynamics of the poles. Substituting (8) in (5) we 
derive the following ordinary differential equations for the 
positions of the poles: 



~dt 



2N 

■ E 



cot 



Zj z k 



i—sign[Im(zj )] 



(10) 



We note that in (8) we have 2N poles which are arranged 
in pairs. In the second sum in (8) each pair of poles 
contributed one term. We can write the pole dynamics in 



A. Qualitative properties of the stationary solution 

The steady-state solution of the flame front propagat- 
ing in channels of width 2tt was presented in Ref. [9]. 
Using these results we can immediately translate the dis- 
cussion to a channel of width L. The main results are 
summarized as follows: 

1. There is only one stable stationary solution which is 
geometrically represented by a giant cusp (or equiv- 
alently one finger) and analytically by N(L) poles 
which are aligned on one line parallel to the imag- 
inary axis. The existence of this solution is made 
clearer with the following remarks. 

2. There exists an attraction between the poles along 
the real line. This is obvious from Eq.(ll) in 
which the sign of dxj/dt is always determined by 
sin ). The resulting dynamics merges all the 
x positions of poles whose ^-position remains finite. 

3. The y positions are distinct, and the poles are 
aligned above each others in positions j/j-i < yj < 
yj+i with the maximal being yj\f(L) ■ This can be 
understood from Eq.(12) in which the interaction 
is seen to be repulsive at short ranges, but changes 
sign at longer ranges. 

4. If one adds an additional pole to such a solution, 
this pole (or another) will be pushed to infinity 
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along the imaginary axis. If the system has less 
than N(L) poles it is unstable to the addition of 
poles, and any noise will drive the system towards 
this unique state. The number N(L) is 



N(L) 



2 V v 



(13) 



where ^---J is the integer part. To see this con- 
sider a system with N poles and such that all the 
values of yj satisfy the condition < yj < y m ax ■ 
Add now one additional pole whose coordinates are 
z a = (x a ,y a ) with y a > y max . From the equa- 
tion of motion for y a , (12) we see that the terms in 
the sum are all of the order of unity as is also the 
cot(t/ a ) term. Thus the equation of motion of y a is 
approximately 



dy a 
dt 



2N + 1 



1_ 

L 



(14) 



The fate of this pole depends on the number of 
other poles. If N is too large the pole will run 
to infinity, whereas if N is small the pole will be 
attracted towards the real axis. The condition for 
moving away to infinity is that N > N(L) where 
N(L) is given by (13). On the other hand the y 
coordinate of the poles cannot hit zero. Zero is 
repulsive line, and poles are pushed away from zero 
with infinite velocity. To see this consider a pole 
whose yj approaches zero. For any finite L the 
term coth(t/j) grows unboundedly whereas all the 
other terms in Eq.(12) remain bounded. 

5. The height of the cusp is proportional to L. The 
distribution of positions of the poles along the line 
of constant x was worked out in [9]. 

We will refer to the solution with all these properties as 
the Thual-Frisch-Henon (TFH)-cusp solution. Next we 
turn to the stability analysis of this solution. 



III. LINEAR STABILITY ANALYSIS IN 
CHANNEL GEOMETRY 

In this section we discuss the linear stability of the 
TFH-cusp solution. To this aim we first use Eq.(8) to 
write the steady solution u s (6) in the form: 



N 



2 sin[0 ■ 



cosh[t/j] — cosf^ — x s ] 



(15) 



where x s is the real (common) position of the station- 
ary poles and yj their stationary imaginary position. To 
study the stability of this solution we need to determine 



the actual positions yj . This is done numerically by in- 
tegrating the equations of motion for the poles starting 
from N poles in initial positions and waiting for relax- 
ation. Next one perturbs this solution with a small per- 
turbation <j){6,t): u(9,i) = u s (6) + <j){6,t) . Linearizing 
the dynamics for small <f> results in the equation of motion 



d<i>(e,t) 

dt 



2 ] d e [u,(9)<f>(9,t)] + vdU(0,t)\ 



I(<K6,t)) 



(16) 



A. Fourier decomposition and eigenvalues 

The linear equation can be decomposed in Fourier 
modes according to 



0= J2 ^*(*) £ 



N 



,{0) = -2vi J2 



I k I y j ikO 

e 1 |yj e 



(17) 
(18) 



-OO 7=1 



In these sums the discrete k values run over all the inte- 
gers. Substituting in (16) we get the equations: 



dt 



y^^akn^njt) 



(19) 



where a^n is a infinite matrix whose entries are given by 



a-kk 



(20) 



L L 2 



N 



To solve for the eigenvalues of this matrix we need to 
truncate it at some cutoff vector k* . The choice of k* 
can be based on the linear stability analysis of the flat 
front. The scale k max , cf. (7), is the largest k which is 
still linearly unstable. We must choose k* > k max and 
test the choice by the converegence of the eigenvalues. 
The chosen value of k* in our numerics was Ak max . The 
results for the low order eigenvalues of the matrix a kn 
that were obtained from a converged numerical calcula- 
tion are presented in Fig.l. The eigenvalues are multi- 
plied by L 2 and are plotted as a function of L. We order 
the eigenvalues in decreasing order and denote them as 
Ao < Ai < A2 . . .. The figure offers a number of qualita- 
tive observations: 

1. There exists an obvious Goldstone or translational 
mode u'(6) with eigenvalue Ao = 0. This eigenmode 
stems from the Galilean invariance of the equation 
of motion. 
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2. The eigenvalues oscillate periodically between val- 
ues that are ^-independent in this presentation (in 
which we multiply by L 2 ). In other words, up to 
the oscillatory behavior the eigenvalues depend on 
L like L~ 2 . 

3. The eigenvalues Ai and A2 hit zero periodically. 
The functional dependence in this presentation ap- 
pears almost piecewise linear. 

4. The higher eigenvalues also exhibit similar qualita- 
tive behaviour, but without reaching zero. We note 
that the solution becomes marginally stable for ev- 
ery value of L for which the eigenvalues Ai and A2 
hit zero. The L~ 2 dependnece of the spectrum in- 
dicates that the solution becomes more and more 
sensitive to noise as L increases. 



B. Qualitative understanding using pole-analysis 

The most interesting qualitative aspects are those enu- 
merated above as item 2 and 3. To understand them it is 
useful to return to the pole description, and to focus on 
Eq.(14). This equation describes the dynamics of a sin- 
gle far-away pole. We remarked before that this equation 
shows that for fixed L the stable number of poles is the 
integer part (13). Define now the number a, < a < 1, 
according to 



-i L - 



-i L - 



Using this number we rewrite Eq.(14) as 



dy a 
dt 



2v 
l 2 ' 



(22) 



(23) 



As L increases, a oscillates piecewise linearly and period- 
ically between zero and unity. This shows that a distant 
pole which is added to the giant cusp solution is usu- 
ally repelled to infinity except when a hits zero when the 
system becomes marginally unstable to the addition of a 
new pole. 

To connect this to the linear stability analysis we note 
from Eq.(8) that a single far-away pole solution (i.e with 
y very large) can be written as 



u(6,t) = Ave 



-y(*)« 



i(0 - x{t)) 



(24) 



Suppose that we add to our giant cusp solution a pertur- 
bation of this functional form. ^From Eq.(23) we know 
that y grows linearly in time, and therefore this solution 
decays exponentially in time. The rate of decay is a lin- 
ear eigenvalue of the stability problem, and from Eq.(23) 
we understand both the 1/L 2 dependence and the peri- 
odic marginality. We should note that this way of think- 
ing gives us a significant part of the L dependence of the 



eigenvalues, but not all. The variable a is rising from zero 
to unity periodically, but after reaching unity it hits zero 
instantly. Accordingly, if the highest non zero eigenvalue 
were fully determined by the pole analysis, we would ex- 
pect this eigenvalue to behave as the dotted line shown 
in Fig. 2. The actual highest eigenvalue is shown in con- 
tinuous line. It is clear that the pole analysis gives us a 
great deal of qualitative and quantitative understanding, 
but not all the features agree. 



C. Dynamics near marginality 

The discovery of marginality at isolated values of L 
poses questions regarding the fate of poles that are added 
at very large y's at certain ^-positions. We will argue 
now that when the system becomes marginally stable, a 
new pole can be added to those existing in the giant cusp. 
We remember that these poles have a common 6 position 
that we denote as 6 = 6 C . The fate of a new pole added 
at infinity depends on its 6 position. If the position of the 
new pole is again denoted as y a , and 00 y a y m ax, 
we can see from Eq.(12) that dy a /dt is maximal when 
9 a = 6 C , whereasit is minimal when a — 6 C = ir. This 
follows from the fact that the cosine term has a value +1 
when 6 a = 6 C and a value —1 when a —6 c = ir. For large 
y differences the terms in the sum take on their minimal 
value when the cos term is —1 and their maximal values 
at +1. For infinitely large y a the equation of motion is 
(14) which is independent of 9 a . Since the RHS of this 
equation becomes zero at marginality, we conclude that 
for very large y a which is not infinite dy a /dt changes sign 
from positive to negative when 9 a — 6 C changes from zero 
to 7r. The meaning of this observation is that the most 
unstable points in the system are those points which are 
furthest away from the giant cusp. It is interesting to 
discuss the fate of a pole that is added to the system 
at such a position. From the point of view of the pole 
dynamics 6 = 6 C + 7r is an unstable fixed point for the 
motion along the 6 axis. The attraction to the giant 
cusp exactly vanishes at this point. If we start with a 
pole at a very large y a close to this value of 6 the down- 
fall along the the y coordinates will be faster then the 
lateral motion towards the giant cusp. We expect to see 
therefore the creation of a small cusp at 6 values close 
to 7r that precedes a later stage of motion in which the 
small cusp moves to merge with the giant cusp. Upon 
the approach of the new pole to the giant cusp all the 
existing poles will move up and the furthest pole at y max 
will be kicked off to infinity. We will later explain that 
this type of dynamics occurs in stable systems that are 
driven by noise. The noise generates far away poles (in 
the imaginary direction) that get attracted around 9 = 
O c + ir to create small cusps that run continuously towards 
the giant cusp. 
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D. Nonlinear Stability 

The intuition gained so far can be used to discuss the 
issue of stability of a stable system to larger perturba- 
tions. In other words, we may want to add to the system 
poles at finite values of y and ask about their fate. We 
first show in this subsection that poles whose initial y 
value is below y max ~ log(_L 2 /V 2 ) will be attracted to- 
wards the real axis. The scenario is similar to the one 
described in the last paragraph. 

Suppose that we generate a stable system with a giant 
cusp at 6 C = with poles distributed along the y axis up 
to y m ax ■ We know that the sum of all the forces that act 
on the upper pole is zero. Consider then an additional 
pole inserted in the position (y max , t). It is obvious from 
Eq.(12) that the forces acting on this pole will pull it 
downward. On the other hand if its initial position is 
much above y max the force on it will be repulsive towards 
infinity. We see that this simple argument identifies y max 
as the typical scale for nonlinear instability. 

Next we estimate y max and interpret our result in 
terms of the amplitude of a perturbation of the flame 
front. To this aim we employ the result of ref. [9] re- 
garding the stable distribution of pole positions in a sta- 
ble large system. The parametrization of [9] differs from 
ours; to go from our parametrization in Eq.(5) to theirs 
we need to rescale u by L~ x and t by L. The parameter v 
in their parameterization is v / L in ours. The distribution 
reads 



p(y) = -^m[coth(M/4)] 



(25) 



To estimate the position of the uppermost pole we require 
that 



dyp(y) 



(26) 



Expanding (25) for large y and integrating explicitly the 
result in (26) we end up with the estimate 



Vn 



2 In 



AL 



(27) 



For large L this result is y max ~ ln(-^2-). If we now add 
an additional pole in the position (y max , 0) this is equiv- 
alent to perturbing the solution u(6,i) with a function 
i>e~ ,Jma * sm(9)/L, as can be seen directly from (8). We 
thus conclude that the system is unstable to a perturba- 
tion larger than 



u{9) ~ v 3 $m(Q)/L 2 



(28) 



This indicates a very strong size dependence of the sensi- 
tivity of the giant cusp solution to external perturbations. 
This will be an important ingredient in our discussion of 
noisy systems. 



IV. INITIAL CONDITIONS, POLE 
DECOMPOSITION AND COARSENING 

In this section we show first that any initial conditions 
can be approximated by pole decomposition. Later, we 
show that the dynamics of sufficiently smooth initial data 
can be well understood from the pole decomposition. Fi- 
nally we employ this picture to describe the inverse cas- 
cade of cusps into the giant cusp which is the final steady 
state. By inverse cascade we mean a nonlinear coarsen- 
ing process in which the small scales coalesce in favor 
of larger scales and finally the system staturates at the 
largest available scale [16]. 



A. Pole Expansion: General Comments 

The fundamental question is how many poles are 
needed to describe any given initial condition. The an- 
swer, of course, depends on how smooth are the initial 
conditions. Suppose also that we have an initial function 
u(6,t = 0) that is 27r-periodic and which at time t = 
admits a Fourier representation 



u(6) = J2 Ak s™(k9 + <j>) 



(29) 



k = l 



with Ak > for all k. Suppose that we want to find a 
pole-decomposition representation u p (6) such that 



\ u p(@) ~ u {6)\ < e f° r every 9 



(30) 



where e is a given wanted accuracy. If u(9) is differen- 
tiable we can cut the Fourier expansion at some finite 
k = K knowing that the remainder is smaller than, say, 
e/2. Choose now a large number M and a small number 
A< 1/M and write the pole representation for u p (6) as 



K M-l 



2ksm(k6 + <f>) 



k = l p = 



cosh [k(yk + pA)] — cos(k6 + <f>) 



(31) 



To see that this representation is a particular form of the 
general formula (8) We use the following two identities 



E-kt 
e > 



"k = — 



k = 



2 cosh t 



(32) 



K-l 



sm(x + ky) = sin(a 



k = 



K - 1 , . Ky y 

2/) sin— cosec- . (33) 



^,From these follows a third identity 



K-l 



E 



2 sin (x 



K 



q cosh y — cos (x jp- + <j>) 

2K sin (Kx + <f>) 
cosh Ky — cos (Kx + <f>) 



(34) 
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Note that the LHS of (34) is of the form (8) with K poles 
whose positions are all on the line yj = y and whose xj 
are on the lattice points 2irj/K — <f>. On the other hand 
every term in (31) is of this form. 

Next we use (32) to rewrite (31) in the form 

K M-l co 

u p (6) = E E E 4ke- nk ^ k+p ^ sin (nk6 + n<f>) . 

k = l p = n = l 

(35) 

Exchanging order of summation between n and p we can 
perform the geometric sum on p. Denoting 



K 



M-l 

*y ^ ^—nkpA 
p = 



1 



-MknA 



-knA 



(36) 



we find 



K 



u p (6) 



E ^ kh n,ke~ nkyk sin (nk9 + n<f>) 

k=l n = l 

K oo 

E 4kb n,ke~ nkyk sin (nk9 + n<f>) 

k=l n = 2 
K 

Y^kb ltk e~ kyk sm(k6 + n<f>) . 



(37) 



k = l 



Compare now the second term on the RHS of (37) with 
(29). We can identify 



-ky k 



4kb lik 

The first term can be then bound from above as 

K oo 



E ikb n,ke~ nkyk sin (nk9 + n<j>) 



k=l n = 2 

K oo 



(38) 



(39) 



<EE 4M - 



k=l n = 2 



sin (nkO + n<f>) 



The sine function and the factor (4A') 1 " can be replaced 
by unity and we can bound the RHS of (39) by 



K oo 

EE 

k=l n=2 



Ah, 



K 



b n ,k < y^,A k 

k=l n = l 



Ah 



(40) 



where we have used the fact that b n k < b\ k which follows 
directly from (36). Using now the facts that b\K < b i,k 
for every k < K and that A k is bounded by some finite 
C since it is a Fourier coefficient, we can bound (40) by 
C 2 K/(bi K — C). Since we can select the free parameters 
A and M to make b\ K as large as we want, we can make 
the remainder series smaller in absolute value than e/2. 

The conclusion of this demonstration is that any ini- 
tial condition that can be represented in Fourier series 



can be approximated to a desired accuracy by pole- 
decomposition. The number of needed poles is of the 
order K x M . Of course, the number of poles thus gen- 
erated by the initial conditions may exceed the number 
N(L) found in Eq.(13). In such a case the excess poles 
will move to infinity and will become irrelevant for the 
short time dynamics. Thus a smaller number of poles 
may be needed to describe the state at larger times then 
at t = 0. We need to stress at this point that the pole 
decomposition is over complete; for example, if there is 
exactly one pole at t = and we use the above tech- 
nique to reach a pole decomposition we would get a large 
number of poles in our representation. 



B. The initial stages of the front evolution: the 
exponential stage and the inverse cascade 

In this section we employ the connection between 
Fourier expansion and pole decomposition to understand 
the initial exponential stage of the evolution of the flame 
front with small initial data u(6,t = 0). Next we employ 
our knowledge of the pole interactions to explain the slow 
dynamics of coarsening into the steady state solution. 

Suppose that initially the expansion (29) is available 
with all the coefficients A k <C 1. We know from the linear 
instability of the flat flame front that each Fourier com- 
ponent changes exponentially in time according to the 
linear spectrum (6). The components with wave vector 
larger than (7) decrease, whereas those with lower wave 
vectors increase. The fastest growing mode is k c = L/2v. 
In the linear stage of growth this mode will dominate the 
shape of the flame front, i.e. 



u(6,t) « At^"^ sm(k c e) 



(41) 



Using Eq.(34) for a large value of y (which is equivalent 
to small A kc ) we see that to the order of 0(A\ ) (41) can 
be represented as a sum over L/2v poles arranged peri- 
odically along the 9 axis. Other unstable modes will con- 
tribute similar arrays of poles but at much higher values 
of y, since their amplitude is exponentially smaller. We 
note that even when all the unstable modes are present, 
the number of poles which they introduce is finite for fi- 
nite L, since there are only L/v unstable modes. Count- 
ing the number of poles that each mode introduces we 

get a total number of (^L/v^ poles. The number of poles 

L/2v which are associated with the most unstable mode 
is precisely the number allowed in the stable stationary 
solution, cf.(13). When the poles approach the real axis 
and cusps begin to develop, the linear analysis no longer 
holds, but the pole description does. 

We now describe the qualitative scenario for the estab- 
lishment of the steady state. Firstly, we understand that 
all the poles that belong to less unstable modes will be 
pushed towards infinity. To see this think of the system 
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at this stage as an array of uncoupled systems with a 
scale of the order of unity. Each such system will have 
a cahracteristic value of y. As we discussed before poles 
that are further away along the y axis will be pushed 
to infinity. Therefore system will remain with the L/2v 
poles of the most unstable mode. The net effect of the 
poles belonging to the (nonlinearly) stable modes is to 
destroy the otherwise perfect periodicity of the poles of 
the unstable mode. 

Once the poles of the stable modes get sufficiently far 
from the real axis, the dynamics of the remaining poles 
will begin to develop according to the interactions that 
are directed along the real axis. These interactions are 
much weaker and the resulting dynamics occur on much 
longer time scales. The qualitative picture is of an in- 
verse cascade of merging the 9 positions of the poles. We 
note that the system has a set of unstable fixed points 
which are 'cellular solutions' described by a periodic ar- 
rangement of poles along the real axis with a frequency 
k. These fixed points are not stable to perturbations and 
they collapse under perturbations with a charachtertic 
time scale (that depends on k) to the next unstable fixed 
point at k' = k/2. This process then goes on indefinitely 
until k ~ l/L i.e. we reach the giant cusp, the steady- 
state stable solution [16]. 

This scenario is seen very clearly in the numerical sim- 
ulations. In Fig. 3 we show the time evolution of the 
flame front starting from small white-noise initial condi- 
tions. The top curve pertains to the earliest time in this 
picture, just after the fast exponential growth, and one 
sees clearly the periodic array of cusps that form. The 
succeeding images show the progress of the flame front in 
time, and one observes the development of larger scales 
with deeper cusps that represent the partial coalescence 
of poles onto the same 6 positions. In Fig. 4 we show the 
width of this front as a function of time. One recognizes 
the exponential stage of growth in which the L/2v poles 
approach the 9 axis, and then a clear cross-over to much 
slower dynamics in which the effective scale in the system 
grows with a slower rate. 

The slow dynamics stage can be understood qualita- 
tively using the previous understanding of the cascade as 
follows: if the initial number of poles belonging to the 
unstable mode is L/2v, the initial effective linear scale 
is 2v . Thus the first step of the inverse cascade will be 
completed in a time scale of the order of 2v . At this point 
the effective linear scale doubles to Av, and the second 
step will be completed after such a time scale. We want 
to know what is the typical length scale l t seen in the 
system at time t. The typical width of the system at this 
stage will be proporional to this scale. 

Denote the number of cascade steps that took place 
until this scale is achieved by s;. The total time elapsed, 
t{lt) is the sum 



*(/ t )~X>'- (42) 

8 = 1 

The geometric sum is dominated by the largest term 
and we therefore estimate t(l t ) ~ It- We conclude that 
the scale and the width are linear in the time elapsed 
from the initial conditions (7 t ~ t'' , ( = 1). In noise- 
less simulations we find (see Fig. 4) a value of £ which is 
C « 0.95 ± 0.1. 

C. Inverse cascade in the presence of noise 

An interesting consequence of the discussion in the last 
section is that the inverse cascade process is an effective 
"clock" that measures the typical time scales in this sys- 
tem. For future purposes we need to know the typical 
time scales when the dynamics is perturbed by random 
noise. To this aim we ran simulations following the in- 
verse cascade in the presence of external noise. The main 
result that will be used in later arguments is that now 
the appearance of a typical scale l t occurs not after time 
t, but rather according to 

l t ~t c , ( ~ 1.2 ±0.1 . (43) 

The numerical confirmation of this law is exhibited in 
Fig. 5 . We also find that the front velocity in this case 
increases with time according to 

v-f 1 , 7 « 0.48 ±0.05 . (44) 

This result will be related to the acceleration of the flame 
front in noisy simulation, as will be seen in the next Sec- 
tions. The result (43) will be helpful in Section V C in 
estimating the values of the scaling exponents. 

V. ACCELERATION OF THE FLAME FRONT, 
POLE DYNAMICS AND NOISE 

A major motivation of this Section is the observation 
that in radial geometry the same equation of motion 
shows an acceleration of the flame front. The aim of 
this section is to argue that this phenomenon is caused 
by the noisy generation of new poles. Moreover, it is our 
contention that a great deal can be learned about the 
acceleration in radial geometry by considering the effect 
of noise in channel growth. In Ref. [9] it was shown that 
any initial condition which is represented in poles goes to 
a unique stationary state which is the giant cusp which 
propagates with a constant velocity v = 1/2 up to small 
l/L corrections. In light of our discussion of the last sec- 
tion we expect that any smooth enough initial condition 
will go to the same stationary state. Thus if there is no 
noise in the dynamics of a finite channel, no acceleration 
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of the flame front is possible. What happens if we add 
noise to the system? 

For concreteness we introduce an additive white-noise 
term rj(9,t)to the equation of motion (5) where 



(45) 



and the Fourier amplitudes rjk are correlated according 
to 



< m (t) 11 l,(t') >=f6(k-k')6(t-t') 



(46) 



We will first examine the result of numerical simulations 
of noise-driven dynamics, and later return to the theo- 
retical analysis. 



A. Noisy Simulations 

Previous numerical investigations [6,15] did not intro- 
duce noise in a controlled fashion. We will argue later 
that some of the phenomena encountered in these sim- 
ulations can be ascribed to the (uncontrolled) numeri- 
cal noise. We performed numerical simulations of Eq.(5) 
with an additive white noise by choosing rjk for every k 
from a flat distribution in the interval [— s/Tf , We 
examined the average steady state velocity of the front 
as a function of L for fixed / and as a function of / for 
fixed L. We found the interesting phenomena that are 
summarized here: 

1. When the noise level / is fixed the average veloc- 
ity v increases with L, see Fig. 6. There are two 
different regimes of this behavior. For sufficiently 
small values of / and L we observe a power law 
dependence of v on L: 



L», ^k, 0.44 ±0.03 



(47) 



For large values of / and L the dependence of v on L 
is much stronger, \i > 1.5, but we cannot estimate 
it well because of lack of dynamical range. For 
a given value of / the transition between the two 
regimes appears upon increasing L, and we denote 
the critical value as L c (f), see Fig. 7. The data 
indicates that 



Lc(f) ~ r 



(48) 



with a « 0.5 ± 0.1. 



2. When the system size L is fixed the average velocity 
depends on / as 



(49) 



For sufficiently small value of / this dependence is 
very weak, and £ ?a 0.02, see Fig. 8. For large values 
of / the dependence is much stronger, £ = 1 ± 0.1. 



3. The locus separating weaker dependencies of v on 
L and / (regime I) from stronger dependencies 
(regime II) is shown in Fig. 7. The giant cusp re- 
mains recognizable in both regimes. In regime I 
one observes a series of small cusps superposed on 
the structure of the giant cusp. In Fig. 9 we show 
a set of interfaces with a growing size and with the 
same noise. One can observe a strong increase in 
the number of cusps and the complexity of their ar- 
rangements as the system size increases. In regime 
II there are strong fluctuations in the gradient field 
u(9,t), see Fig. 10. 

4. Measurements of the width of the front indicate 
that it has a a very weak dependence on / for fixed 
L in regime I. 



B. Theoretical Discussion of the Effect of Noise 

1. The Threshold of Instability to Added Noise 

First we present the theoretical arguments that explain 
the sensitivity of the giant cusp solution to the effect of 
added noise. This sensitivity increases dramatically with 
increasing the system size L. To see this we use again 
the relationship between the linear stability analysis and 
the pole dynamics. 

Our additive noise introduces perturbations with all 
fc-vectors. We showed previously that the most unsta- 
ble mode is the k = 1 component Aism(9). Thus the 
most effective noisy perturbation is r\\ sin(#) which can 
potentially lead to a growth of the most unstable mode. 
Whether or not this mode will grow depends on the am- 
plitude of the noise. To see this clearly we return to the 
pole description. For small values of the amplitude A\ we 
represent A\ sin(#) as a single pole solution of the func- 
tional form (y I ' L)e~ y sm9. The y position is determined 
from y = — log \Ai\L/v, and the ^-position is 9 = 7r for 
positive A\ and 9 = for negative Ai . From the analysis 
of Section III we know that for very small Ai the fate 
of the pole is to be pushed to infinity, independently of 
its 9 position; the dynamics is symmetric in Ai — ► — Ai 
when y is large enough. On the other hand when the 
value of Ai increases the symmetry is broken and the 9 
position and the sign of Ai become very important. If 
Ai > there is a threshold value of y below which the 
pole is attracted down. On the other hand if Ai < 0, 
and 9 = the repulsion from the poles of the giant cusp 
grows with decreasing y. We thus understand that qual- 
itatively speaking the dynamics of Ai is characterized by 
an asymmetric "potential" according to 



A 1 



dAi 



ViA,) = \A\ 



aA\ 



(50) 
(51) 
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^From the linear stability analysis we know that A ps 
v I ' L 2 , cf. Eq.(14). We know further that the threshold 
for nonlinear instability is at A\ ps v 3 j ' L 2 , cf. Eq(28). 
This determines that value of the coefficient a ps 2/3v 2 . 
The magnitude of the "potential" at the maximum is 



V(A max ) « v r /L 6 



(52) 



The effect of the noise on the development of the mode 
Ai sin 6 can be understood from the following stochastic 
equation 



A! 



dA 1 



(53) 



It is well known [17] that for such dynamics the rate of 
escape R over the "potential" barrier for small noise is 
proportional to 



R ■ 



v 

L 2 



(54) 



The conclusion is that any arbitrarily tiny noise becomes 
effective when the system size increase and when v de- 
creases. If we drive the system with noise of amplitude 
/ the system can always be sensitive to this noise when 
its size exceeds a critical value L c that is determined by 
/ ~ L & c /v 7 . For L > L c the noise will introduce new 
poles into the system. Even numerical noise in simula- 
tions involving large size systems may have a macroscopic 
influence. 

The appearance of new poles must increase the veloc- 
ity of the front. The velocity is proportional to the mean 
of (u/L) 2 . New poles distort the giant cusp by additional 
smaller cusps on the wings of the giant cusp, increasing 
u 2 . Upon increasing the noise amplitude more and more 
smaller cusps appear in the front, and inevitably the ve- 
locity increases. This phenomenon is discussed quantita- 
tively in Section V. 



2. The Noisy Steady State and its Collapse with Large Noise 
and System Size 

In this subsection we discuss the response of the giant 
cusp solution to noise levels that are able to introduce a 
large number of excess poles in addition to those exist- 
ing in the giant cusp. We will denote the excess number 
of poles by 6N . The first question that we address is 
how difficult is it to insert yet an additional pole when 
there is already a given excess 6N . To this aim we esti- 
mate the effective potential Vsn(Ai) which is similar to 
(51) but is taking into account the existence of an excess 
number of poles. A basic approximation that we employ 
is that the fundamental form of the giant cusp solution 
is not seriously modified by the existence of an excess 
number of poles. Of course this approximation breaks 



down quantitatively already with one excess pole. Qual- 
itatively however it holds well until the excess number of 
poles is of the order of the original number N(L) of the 
giant cusp solution. Another approximation is that the 
rest of the linear modes play no role in this cse. At this 
point we limit the discussion therefore to the situation 
6N < N(L). 

To estimate the parameter A in the effective potential 
we consider the dynamics of one pole whose y position 
y a is far above y ma x- According to Eq.(14) the dynamics 
reads 



dju _ v(N(L) + 8N) 1 
dt ~ L 2 L 



(55) 



Since the N(L) term cancels against the L~ x term (cf. 
Sec. II A), we remain with a repulsive term that in the 
effective potential translates to 



A 



v8N 



(56) 



Next we estimate the value of the potential at the break- 
even point between attraction and repulsion. In the last 
subsection we saw that a foreign pole has to be inserted 
below y m ax m order to be attracted towards the real axis. 
Now we need to push the new pole below the position 
of the existing pole whose index is N(L) — 8N. This 
position is estimated as in Sec III C by employing the 
TFH distribution function (25). We find 



y 6N ps 2 In 



4L 



Tt 2 vbN 



(57) 



As before, this implies a threshold value of the amplitude 
of single pole solution A max sin 9 which is obtained from 
equating A max = ve y6N . We thus find in the present case 
Amax ~ v 3 (6N) 2 / L 2 . Using again a cubic representation 
for the effective potential we find a = 2/(36N) and 



V (A max ) 



lv 7 {6Nf 
3 U~ 



(58) 



Repeating the calculation of the escape rate over the po- 
tential barrier we find in the present case 



R 



L 2 



■ exp 



- v r (6N) 5 /fL 6 



(59) 



For a given noise amplitude / there is always a value 
of L and v for which the escape rate is of 0(1) as long 
as 6N is not too large. When 6N increases the escape 
rate decreases, and eventually no additional poles can 
creep into the system. The typical number 6N for fixed 
values of the parameters is estimated from equating the 
argument in the exponent to unity: 



(60) 



The most important consequence of this relation is that 
6N increases with L faster than N(L). Accordingly we 
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expect a breakdown of this picture and of the weak noise 
behavior when 6N ps N(L), which is occurring when L 
reaches a critical value L c (f) where 

Mi) ~ r 1 • (6i) 

As before, this prediction is not in quantitative agree- 
ment with (48) for the reasons discussed above. Never- 
theless we believe that it provides a reasonable qualita- 
tive understanding for the origin of the power law behav- 
ior that is found numerically. 

C. The acceleration of the flame front due to noise 

In this section we estimate the scaling exponents that 
characterize the velocity of the flame front as a function of 
the system size. Our arguments in this section are even 
less solid than the previous ones, but nevertheless we 
believe that we succeed to capture some of the essential 
qualitative physics that underlies the interaction between 
noise and instability and which results in the acceleration 
of the flame front. 

To estimate the velocity of the flame front we need 
to write down an equation for the mean of < dh/dt > 
given an arbitrary number N of poles in the system. This 
equation follows directly from (4): 

/ dh\ 1 1 f 2w , ,„ 

After substitution of (8) in (62) we get, using (11) and 
(12) 




The estimates of the second and third terms in this equa- 
tion are straightforward. Writing N = N(L) + 8N(L) 
and remembering that N(L) = v / L we find that these 
terms contribute —vhN(L)/L ~ —L 1 / 5 . The first term 
will contribute only when the current of poles is asym- 
metric. Since noise introduces poles at a finite value of 
y, whereas the rejected poles stream towards infinity, we 
have an asymmetry that contributes to the velocity of the 
front. To estimate the first term we note that in noiseless 
conditions the typical time in this system scales like L, 
whereas in the presence of noise (cf. Eq.(43)) like L~ . 
Accordingly, the typical flux of poles can be estimated 
as 8N(L)/L~ . Thus the current ^2 Vk has a stronger de- 
pendence on L, i.e _L 6 / 5-1 /^ . Taking the numerical value 
of C = 1.2 we conclude that Eq.(63) predicts a scaling 
law (47) with /j, = 0.37, in reasonable agreement with 
the numerics. 

We should stress at this point that the argument is 
not complete. We measured the exponent relating the 



velocity to the amplitude of external noise, cf. Eq.(49) 
and found that £ is considerably smaller than the value 
1/5 which is predicted by the previous argument. This 
indicates that the typical time scale also has an explicit 
/ dependence which is physically plausible but it was not 
measured by us. 



VI. SUMMARY AND CONCLUSIONS 

The main two messages of this paper are: (i) There is 
an important interaction between the instability of de- 
veloping fronts and random noise; (ii) This interaction 
and its implications can be understood qualitatively and 
sometime quantitatively using the description in terms of 
complex poles. 

The pole description is natural in this context firstly 
because it provides an exact (and effective) representa- 
tion of the steady state without noise. Once one succeeds 
to describe also the perturbations about this steady state 
in terms of poles, one achieves a particularly transparent 
language for the study of the interplay between noise and 
instability. This language also allowed us to describe in 
qualitative and semi-quantitative terms the inverse cas- 
cade process of increasing typical lengths when the sys- 
tem relaxes to the steady state from small, random initial 
conditions. 

The main conceptual steps in this paper were as fol- 
lows: firstly one realizes that the steady state solution, 
which is characterized by N(L) poles aligned along the 
imaginary axis is marginally stable against noise in a pe- 
riodic array of L values. For all values of L the steady 
state is nonlinearly unstable against noise. The main and 
foremost effect of noise of a given amplitude / is to intro- 
duce an excess number of poles 8N(L, /) into the system. 
The existence of this excess number of poles is responsi- 
ble for the additional wrinkling of the flame front on top 
of the giant cusp, and for the observed acceleration of the 
flame front. By considering the noisy appearance of new 
poles we could rationalize the observed scaling laws as a 
function of the noise amplitude and the system size. 

Theoretically we therefore concentrate on estimating 
8N(L, /). The measurements are not testing our theoret- 
ical consideration directly, but rather test the dependence 
of the velocity on L and /. The only direct test is the 
critical line shown in Fig. 7, and there is a significant dif- 
ference between our predictions and the measurements. 
We understand that some of our consideration are only 
qualitative. For example, we estimated 8N(L,f) by as- 
suming that the giant cusp solution is not seriously per- 
turbed. On the other hand we found a flux of poles going 
to infinity due to the introduction of poles at finite values 
of y by the noise. The existence of poles spread between 
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Umax and infinity is a significant perturbation of the gi- 
ant cusp solution. Thus also the comparison between the 
various scaling exponents measured and predicted must 
be done with caution; we cannot guarantee that those 
cases in which our prediction hit close to the measure- 
ment mean that the theory is quantitative. Nevertheless 
we believe that our consideration extract the essential 
ingredients of a correct theory. 

The "phase diagram" as a function of L and / in this 
system consists of three regimes. In the first one, dis- 
cussed in Section V B 1 , the noise is too small to have 
any effect on the giant cusp solution. In the second the 
noise introduces excess poles that serve to decorate the 
giant cusp with side cusps. In this regime we find scaling 
laws for the velocity as a function of L and / and we are 
reasonably successful in understanding the scaling expo- 
nents. In the third regime the noise is large enough to 
create small scale structures that are not neatly under- 
stood in terms of individual poles. It appears from our 
numerics that in this regime the roughening of the flame 
front gains a contribution from the the small scale struc- 
ture in a way that is reminiscent of stable, noise driven 
growth models like the Kardar-Parisi-Zhang model. 

One of our main motivations in this research was to 
understand the phenomena observed in radial geometry 
with expanding flame fronts. A full analysis of this prob- 
lem cannot be presented here. We note however that 
many of the insights offered above translate immediately 
to that problem. Indeed, in radial geometry the flame 
front accelerates and cusps multiply and form a hierar- 
chic structure as time progresses. Since the radius (and 
the typical scale) increase in this system all the time, new 
poles will be added to the system even by a vanishingly 
small noise. The marginal stability found above holds 
also in this case, and the system will allow the intro- 
duction of excess poles as a result of noise. The results 
discussed in Ref. [7] can be combined with the present in- 
sights to provide a theory of radial growth. This theory 
will be offered in a forthcoming publication. 

Finally, the success of this approach in the case of flame 
propagation raises hope that Laplacian growth patterns 
may be dealt with using similar ideas. A problem of 
immediate interest is Laplacian growth in channels, in 
which a finger steady-state solution is known to exist. It 
is documented that the stability of such a finger solution 
to noise decreases rapidly with increasing the channel 
width. In addition, it is understood that noise brings 
about additional geometric features on top of the finger. 
There are enough similarities here to indicate that a care- 
ful analysis of the analytic theory may shed as much light 
on that problem as on the present one. 
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Figures Legends 

Fig.l: The first 10 highest eigenvalues of the stability 
matrix, multiplied by the square of the system size L 2 
vs. the system size L. Note that all the eigennvalues os- 
cillate around fixed values in this presentation, and that 
the highest two eigenvalues hit zero periodically. 

Fig. 2: Comparison of the numerically determined highest 
4 eigenvalues of the stability matrix with the prediction 
of the pole analysis. The latter provides a qualitative un- 
derstanding of the stability, and appears to overlap with 
the highest eigenvector over half of the range, and with 
the fourth eigenvalue over the other half. 

Fig. 3: The inverse cascade process of coarsening that oc- 
curs after preparing the system with random, small initial 
conditions. One sees that at successive times the typical 
scale increases until the giant cusp forms, and attracts all 
the other side-poles. The effect of the existing numerical 
additive noise is to introduce poles that appear as side 
cusps that are continuously attracted to the giant cusp. 
This effect is obvious to the eye only after the typical 
scale is sufficiently large, as is seen in the last time (see 
text for further details). 

Fig. 4: log-log plots of the front velocity (lower curve) 
and width (upper curve) as a function of time in the in- 
verse cascade process seen in Fig. 3 in asystem of size 2000 
and v = 1. Both quantities exhibit an initial exponential 
growth that turns to a power law growth (after t ps 30). 
The velocity is constant after this time, and the width 
increases like . Note that at the earliest time there is a 
slight decrease in the velocity; this is due to the decay of 
linearly stable modes that exist in random initial condi- 
tions. 

Fig. 5: The same as Fig. 4 but with additive random noise- 
for a system of size 1000, v = 0.1 and / = 10 -8 . The 
velocity does not saturate now, and the exponent £ char- 
acterizing the increase of the width with time changes to 
C = 1.2 ± 0.1. The velocity increases in time like t 1 with 
7 ps 0.48 ± 0.04. 

Fig. 6: log-log plot of the velocity as a function of system 
size for different values of the noise amplitude, / = 2.5 x 
10- 13 ,2.5xl0" 9 ,2.5xl0- 7 , 10- 6 ,4xl0- 6 ,2.5xl0" 5 and 
v = 0.1. For low values of / we observe a power law be- 
havior, v ~ , \i ps 0.44 ± 0.03. For larger values of 
/ there is a cross over to a stronger dependence on the 
system size, see text for details. 

Fig. 7: The critical value of L as a function of / for which 
the transition from weaker to stronger dependence of v 
on L is observed, see Fig. 6. In regime I the noise is rela- 
tively small and the cusp picture is qualitatively correct. 
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In regime II the noise is too large to leave the cusp pic- 
ture intact, and see text and Figures 9 and 10 for further 
details. The data indicates a power law as is explained 
in the text. 

Fig. 8: log-log plot of the velocity vs the noise amplitude 
for different system sizes L = 5, 10, 40, 160, v = 0.1. 

Fig. 9: A typical flame fronts in regime I of Fig. 7 where 
the system is sufficiently small not to be terrible affected 
by the noise. The effect of noise in this regime is to add 
additional small cusps to the giant cusp. In figures a-d 
we present fronts for growing system sizes from 10 to 80, 
v = 0.1 and / = 2.5 x 10 -9 . One can observe that when 
the system size grows there are more cusps with a more 
complex structure. 
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Fig. 10: A typical flame front in regime II of Fig. 7. The 
system size is 160, and the noise amplitude is 2.5 x 10 -5 
. This is sufficient to cause a qualitative change in the 
appearance of the flame front: the noise introduces sig- 
nificant levels of small scales structure in addition to the 
cusps. 
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